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Resumen 


Las inundaciones están asociadas con diversos factores como: 1) desbordamiento 
de ríos; 2) inundaciones súbitas, 3) mareas altas asociadas con huracanes, y 
4) rompimiento de estructuras de control. Desafortunadamente, los registros 
históricos de inundación son prácticamente inexistentes. Por lo tanto, los análisis 
de riesgos generalmente se basan en métodos indirectos, como la modelación 
físico-matemática. Se presenta el desarrollo de un Mapa Nacional de un Índice 
de Inundación que permita identificar zonas propensas a inundación. El índice 
propuesto está basado en el Índice Topográfico desarrollado por Beven y Kirkby 
(1979), pero adicionalmente considera factores edafológicos, hidrológicos y 
climatológicos. El proceso se detalla para la Región Hidrológica Lerma-Santiago. 
Adicionalmente, los resultados son comparados con la simulación correspondiente 
de un modelo de tránsito de avenidas detallado de la planicie de inundación de 
Tabasco, México. Los resultados indican que el índice captura razonablemente 
bien las zonas perenes y efímeramente inundadas. Se incluye el Mapa Nacional 
Mexicano del Índice de Inundación obtenido con base en la metodología 
propuesta. 
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Santiago. 


Introducción 


Las inundaciones están asociadas con diversos 
factores como: 1) desbordamiento de ríos; 
2) mareas altas asociadas con huracanes, y 
3) rompimiento de estructuras de control. 
Desafortunadamente, los registros históricos de 
inundación son prácticamente inexistentes. Por 
lo tanto, los análisis de riesgo generalmente se 
basan en métodos indirectos, como la modelación 
físico-matemática. En dicho contexto, la aplicación 
de modelos generalmente tiene como objetivo 
documentar la magnitud y probabilidad de 
daños por inundación, a partir de simulaciones 
basadas en registros históricos de lluvia y caudal. 
Se propone entonces el desarrollo de un Mapa 
Nacional de un Índice de Inundación que permita 
identificar las zonas propensas a inundación con 
base en análisis simplificados y acordes con la 
información disponible en el ámbito nacional. 


El índice propuesto se basa en el Índice 
Topográfico propuesto por Beven y Kirkby 
(1979). El Índice Topográfico es definido como 
el cociente entre acumulación de flujo (área de 
drenaje parcial “aguas arriba” para un punto 
en particular) y la tangente de la pendiente. Su 
aplicación principal consiste en la identifica- 
ción de humedales (wetlands), definidos como 
zonas perennes O efímeramente saturadas 
o inundadas. Debido a que la distribución 
probabilística del Índice Topográfico es 
bimodal para una cuenca, y los valores altos 
del índice corresponden a regiones propensas 
a inundación, el índice constituye un elemento 
poderoso y simple para la determinación de 
dichas regiones (Kirkby et al., 1995; Rohde y 
Seibert, 1999; Merot et al., 1993; entre otros). 
La literatura documenta variaciones del Índice 
Topográfico. Las diferencias radican en la 
implementación de información adicional, co- 
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motipo y uso de suelo, precipitación, humedad 
de suelo, evapotranspiración, etcétera. 

A pesar de que el término “Acumulación de 
Flujo” sugiere que se consideran los volúmenes 
de aguainvolucrados, estono sucedeen realidad, 
pues se trata de una variable derivada a partir 
únicamente de la topografía. Debido a que las 
características climatológicas y de tipo y uso de 
suelo del país son diversas, aunque dos cuencas 
tuvieran características topográficas similares, 
su riesgo de inundación podría ser muy 
distinto. De esa manera, el Índice Topográfico 
resulta insuficiente para caracterizar el riesgo de 
inundación en el país. Por lo cual, se propone un 
índice modificado que permita capturar dichas 
diferencias. 


Datos 


Para la estimación del índice de inundación se 
empleó la siguiente información: 


+ Topografía (continuo de elevación 
digital (CEM) de INEGI con pixel de 
aproximadamente 30 m). 

+ Mapas de vegetación y uso de suelo (INEGI, 
2004, 1:250,000). 

e Mapa de 
1:250,000). 


edafología (INEGI, 2006, 


Adicionalmente, para solventar la falta de 
homogeneidad en cobertura espacial y temporal 
de la base de datos climatológica nacional, 
Agroasemex, S.A., ha desarrollado Arreglos 
Regulares (mallas) de precipitación y otras 
variables climatológicas a partir del método 
de Cressman (1959). Las mallas ofrecen una 
base de datos climatológica completa espacial y 
temporalmente en todo el país para el periodo 
1979-2008. El método de Cressman consiste en la 
corrección de un campo preliminar (estimación 
indirecta de la variable de interés) a partir 
de observaciones. El campo preliminar está 
constituido por los datos del North American 
Regional Reanalysis (NARR) (Mesinger et 
al., 2006) elaborado por el National Center 
for Environmental Prediction del National 


Oceanic and Atmospheric Administration. Las 
observaciones de precipitación son tomadas 
de CLICOM (Clima Computarizado) y GASIR 
(Gerencia de Aguas Superficiales e Ingeniería 
de Ríos) de la Comisión Nacional del Agua 
(Conagua). Las mallas cuentan con una 
resolución temporal diaria y pixeles de 0.2” (-22 


km). 
Metodología 


Se propone un índice definido como el cociente 
entre acumulación de escurrimiento y la 
tangente de la pendiente. El escurrimiento es 
estimado a partir del tipo y uso de suelo, y la 
precipitación observada según el método de 
número de curva. De manera que es posible 
estimar la cantidad total de agua que escurre 
por un punto dado, así como su probabilidad 
de ocurrencia. El análisis se realiza por cada 
una de las cuencas, por lo cual, el primer 
paso consiste en la extracción del CEM de la 
región hidrológica de interés. Los sumideros 
son rellenados y los picos son removidos 
para deponer imperfecciones en los datos. 
Posteriormente se estima la dirección de flujo 
de cada celda a la de su vecina más cercana. 
Además, se calcula la pendiente y su tangente. 
Por otra parte, el escurrimiento se calcula 
mediante el método de curva numérica. 
Debido a que las clasificaciones de tipo de 
suelo de INEGI no son compatibles con las 
tablas de correspondencia de número de curva 
tradicionales, desarrolladas en Estados Unidos 
de Norteamérica, se hizo una reclasificación 
que permitiera la aplicación de las tablas de 
correspondencia. Los valores del número de 
curva bajo condiciones normales de humedad 
de suelo fueron pasados a una retícula para 
cubrir una a una las celdas del CEM. 

Una vez estimado el número de curva 
correspondiente para cada pixel del CEM, se 
procedió a estimar el escurrimiento de acuerdo 
con el método de curva numérica: 


(P-0.25) 


L= 085 
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Donde: 


P = precipitación acumulada (mm). 

S = parámetro deretención máxima de humedad 
del suelo (mm), que depende del número de 
curva (CN), según la siguiente ecuación: 


S= aci 254 
CN 

Debido a que se trata de estimar zonas 
propensas a inundación por causa de eventos 
extremos de precipitación, los valores de 
precipitación empleados para cada región 
corresponden al percentil 95 de cada celda del 
arreglo regular de precipitación de Agroasemex, 
S.A. 

A continuación se estima la acumulación de 
escurrimiento. La diferencia entre acumulación 
de flujo y acumulación de escurrimiento radica 
en que la primera consiste en contar el número 
de celdas que contribuyen a cada celda, mientras 
que en la segunda se cuentan los milímetros de 
escurrimiento asociados con cada celda. 


Por último se estima el índice de inundación 
con base en la siguiente expresión: 


me In ae e] 
tan s 
Donde: 
Il =índice de inundación. 


Ac Q =acumulación de escurrimiento. 
S = pendiente. 


Resultados 


La metodología requiere ser aplicada para 
cuencas completas. Por ello, y por razones de 
simplicidad, eficiencia y presentación, el índice 
de inundación se ha estimado individualmente 
para cada una de las 37 regiones hidrológicas 
que conforman el territorio mexicano (ver figura 
1 y cuadro 1). 

A continuación se analiza la metodología 
y sus resultados en la Región Hidrológica 
Número 12, denominada “Lerma-Santiago”. 


Regiones 


este (Ensenada) 
Oeste (Vizcaino) 


) 


Figura 1. Mapa de las 37 regiones hidrológicas de la república mexicana. 
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Cuadro 1. Regiones hidrológicas de la república mexicana. 


Núm. Nombre de la región Núm. Nombre de la región 
L Baja California Noroeste (Ensenada) 20 Costa Chica y Río Verde 
2 Baja California Centro-Oeste (El Vizcaíno) 21 Costa de Oaxaca (Puerto Ángel) 
3 Baja California Sur-Oeste (Magdalena) 22 Río Tehuantepec 
4 Baja California Noroeste (Laguna Salada) 28) Costa de Chiapas 
5 Baja California Centro-Este (Sta. Rosalía) 24 Río Bravo-Río Conchos 
6 Baja California Sur-Este (La Paz) 25 Ríos San Fernando y Soto La Marina 
7 Río Colorado 26 Río Pánuco 
8 Sonora Norte 27) Ríos Tuxpan y Nautla 
9 Sonora Sur 28 Río Papaloapan 
10 Sinaloa 2 Río Coatzacoalcos 
11 Ríos Presidio y San Pedro 30 Ríos Grijalva-Usumacinta 
(1 Ríos Lerma y Santiago sil Yucatán Oeste (Campeche) 
13 Río Huicicila 32 Yucatán Norte (Yucatán) 
14 Río Ameca 33 Yucatán Este (Quintana Roo) 
15 nd 34 Cuencas Cerradas del Norte 
(Casas Grandes) 
16 Ríos Amería y Coahuayana 39 Mapimí 
17 Costa de Michoacán 36 Ríos Nazas y Aguanaval 
18 Río Balsas 37 El Salado 
19 Costa Grande 


Como la metodología es la misma para el 
resto de las regiones, se considera innecesario 
realizar el mismo análisis para todas ellas. La 
RH12 se encuentra conformada principalmente 
por los estados de Jalisco, Nayarit, Zacatecas, 
Aguascalientes, Guanajuato y, en una menor 
porción. por Durango, Michoacán, Querétaro 
y el Estado de México. Dentro de los cauces 
más importantes se pueden nombrar el río 
Lerma, el río Grande Santiago y el río Verde. 
El almacenamiento o vaso regulador más 
importante con que cuenta la zona es la laguna 
de Chapala. Esta región tiene una superficie 
aproximada de 133 455 km”. 

La figura 2 muestra el CEM de la Región 
Hidrológica Número 12 con cotas que se 
encuentran entre O y 4 660 msnm. 

El resultado de la determinación del 
escurrimiento puede observarse en la 
figura 3; Su rango de variación se encuentra 
entre O y 130 mm. Es posible notar que 
existen dos tipos de variaciones espaciales, 
resultado de arreglos regulares de diferentes 
resoluciones: 1) Precipitación anual con una 


resolución espacial de aproximadamente 
20 km, y 2) Campo de número de curva con 
una resolución espacial de aproximadamente 
30 metros. 

El mapa correspondiente a la tangente de 
la pendiente, es decir, el denominador del 
índice de inundación, se muestra en la figura 
4. Las zonas más oscuras indican pendientes 
pequeñas, asociadas con sitios susceptibles 
de ser inundables por su poca capacidad de 
escurrimiento por gravedad. Las zonas con 
pendientes grandes son susceptibles a generar 
mayor cantidad de escurrimiento y, por lo 
tanto, son relativamente más propensas a 
erosión hídrica. 

Una vez que se obtienen los parámetros 
necesarios para determinar el índice de 
inundación, se estima el índice de inundación 
de acuerdo con la expresión indicada en 
la metodología. El índice de inundación 
resultante es mostrado en la figura 5. Las 
zonas rojas indican una alta susceptibilidad a 
inundación, mientras que las gris claro indican 


lo opuesto. El índice es continuo, y varía 
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Figura 3. Escurrimiento de la RH12. 
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Figura 5. Mapa de índice de inundación de la RH12. 
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entre aproximadamente -14 y 28. Aunque los 
valores del índice pueden ser empleados para 
determinar el riesgo de inundación relativo 
entre diferentes regiones hidrológicas, el análisis 
puede ser simplificado para una misma cuenca 
mediante una clasificación binaria del índice 
de inundación. La clasificación se basa en el 
hecho de que el índice de inundación tienen 
una distribución probabilística bimodal (figura 
6). Siguiendo la convención encontrada en la 
literatura, la porción superior del histograma 
corresponde a regiones susceptibles de ser 
inundadas, mientras que la porción inferior 
corresponde a zonas libres de esta amenaza. 

La figura 7 muestra el resultado de la 
clasificación binaria, indicando en gris oscuro 
las zonas identificadas con riesgo de inun- 
dación. En este mapa es posible observar que 
las zonas propensas a inundación coinciden con 
regiones de baja pendiente, e incluso es posible 
capturar cuerpos de agua importantes, como 
el lago de Chapala y el cauce de algunos ríos. 
Parecería innecesario identificar el riesgo de 


que dichas zonas sufran exceso de humedad, 
pues su naturaleza perenne lo hace evidente; 
sin embargo, es importante recordar que la 
metodología no incluye información alguna 
de la ubicación de dichos cuerpos de agua. 
Por lo tanto, estos resultados sugieren que el 
índice tiene la capacidad de identificar regiones 
perennemente inundadas, y se espera que 
su validez se extienda a zonas efímeramente 
saturadas o cubiertas de agua. La metodología 
anteriormente descrita es aplicada a todas las 
regiones hidrológicas. 

Por último, con el objetivo de corroborar la 
validez de los resultados por otros medios, se 
estimó la erosión hídrica con base en la Ecua- 
ción Universal de Pérdida de Suelo (artículo 
en proceso de publicación, Montes et al., 2009). 
En la figura 8 se muestran los resultados en 
la Región Hidrológica Número 12. Las zonas 
susceptibles de ser inundables coinciden con 
las zonas que presentan un índice de erosión 
bajo. Dicho resultado es el esperado, pues 
en estas zonas se depositan los sedimentos 
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Figura 6. Histograma en donde se puede apreciar que los datos presentan una distribución bimodal. 
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Figura 7. Zonas propensas a inundación por evento extremo de precipitación (gris oscuro) 
en la Región Hidrológica Número 12. 
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Figura 8. Mapa de erosión hídrica de la RH12. 
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producidos en las áreas con altos índices de 
erosión. Cabe mencionar que con el paso del 
tiempo, la acumulación de sedimentos provoca 
que los cauces pierdan capacidad hidráulica, y 
los lagos o lagunas su capacidad de regulación, 
propiciando que se incremente el índice de 
inundación. 


Validación de resultados 


Debido a que se cuenta con un modelo 
hidrológico e hidráulico para una porción 
de la planicie de inundación de Tabasco, 
las simulaciones de dicho modelo fueron 
comparadas con el índice de inundación en 
la misma región, con el objetivo de estimar su 
confiabilidad. La región modelada (figura 9) 
se encuentra dentro de la Región Hidrológica 
Número 30, denominada Grijalva-Usumacinta, 
que comprende varios estados de la república 
mexicana, como Tabasco, Chiapas y Campeche. 
La región cuenta con los ríos más caudalosos 
del país, entre los que se pueden mencionar 
El Grijalva, Usumacinta, Carrizal, Samaria 
y Mezcalapa. Además, tiene el complejo 
hidroeléctrico más grande del país, compuesto 
por cuatro plantas hidroeléctricas: Chicoasén, 
Angostura, Malpaso y Peñitas, lo que permite 
controlar los escurrimientos hacia aguas abajo 
del cauce del río Mezcalapa. 

El modelo hidrológico-hidráulico de la 
cuenca fue desarrollado con anterioridad 
por Agroasemex, S.A., con el propósito 
de desarrollar un esquema de seguro por 
inundación agropecuario en Tabasco. Se trata 
de un modelo detallado en la cuenca de los 
ríos de La Sierra, que comprende los cauces de 
los ríos Tapijulapa, Pichucalco, Teapa, Tepaté, 
La Sierra y Grijalva, cuyas características 
son, entre otras: seccionamiento transversal 
sobre el cauce de los ríos (280 secciones 
transversales); modelos de elevación digital 
escala 1:5,000 hasta 1:20,000; altura de los 
bordos longitudinales existentes, y geometría 
y ubicación de obras hidráulicas. Asimismo, 
se realizó un estudio hidrológico, donde 
se analizaron factores de simultaneidad 


e hidrogramas asociados con diferentes 
periodos de retorno, identificación de las 
zonas susceptibles de ser inundadas y 
conexiones entre éstas. 

Como en la frontera aguas arriba se cuenta 
con hidrogramas asociados con las estaciones 
hidrométricas, y aguas abajo una curva Q-H 
característica de la sección seleccionada, el 
modelo se calibró en las escalas y estaciones 
ubicadas en la parte media y baja de la 
cuenca, de acuerdo con los niveles emitidos 
por Conagua, y se validó la extensión de la 
mancha de afectación con imágenes de satélite 
proporcionadas por el Cenapred, asociadas 
con el evento ocurrido en 2007. Debido a la 
cantidad de información y nivel de detalle, 
es de esperarse que los resultados sean 
diferentes. Sin embargo, se confía en que la 
comparación permita estimar la efectividad 
del índice propuesto para la identificación de 
zonas propensas a inundación. 

La figura 10 muestra la topografía de la 
zona y la extensión de la afectación asociada 
con el evento del año 2007 en la zona de los 
ríos de La Sierra, según el modelo hidrológico- 
hidráulico de Agroasemex, S.A. Se puede 
apreciar la importancia del relieve, pues las 
zonas altas del MDE no son afectadas. Por 
otra parte, la figura 11 muestra en diferentes 
tonos de gris las regiones propensas a 
según 


propuesta (gris oscuro) y el área simulada 


inundación, la metodología aquí 
con el modelo hidrológico hidráulico (gris 
claro). Ambas regiones coinciden en su mayor 
parte de manera satisfactoria, especialmente 
considerando la gran diferencia que existe en 
la información involucrada y metodologías 
empleadas. Las diferencias en las fronteras 
este y oeste se deben a que las fronteras del 
dominio en el modelo hidráulico son más 
limitadas. 


Mapa Nacional de Indice de Inundación: 
consideraciones 


La figura 12 muestra los resultados de aplicar 
la metodología descrita a escala nacional. 
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Figura 9. Ubicación del dominio del modelo hidrológico-hidráulico de Agroasemex, S.A., 


(recuadro en el estado de Tabasco) e hidrografía de la zona. 


Debido a que la resolución del CEM de 
INEGI es de 30 m, el mapa nacional original 
resulta sobresaturado de información; es 
decir, aparentemente todo el país está en 
riesgo de inundación. Por ello, los resultados 
originales son procesados para obtener un 
mapa apropiado para su impresión en la 
actual revista. El procesamiento consiste en 
estimar el análisis de índice de inundación a 


los datos de topografía del SRT'M, cuya menor 
resolución espacial (-90 m) es más apropiada 
para impresión. 

Por su menor resolución, las características 
que refleja el mapa son cuerpos de agua 
mayores, como lagos y planicies de inundación. 
Algunas de las características que son filtradas 
con el aumento de escala son ríos relativamente 
delgados. Sin embargo, el hecho de que sea 
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Figura 11. Mapa de zonas propensas a inundación según el 


índice de inundación (gris oscuro) y mancha de inundación 


según el modelo hidráulico detallado de los ríos de La 


Figura 10. Mancha de inundación asociada con el evento Sierra (gris claro). 


de 2007, según el modelo hidrológico-hidráulico 


de Agroasemex, S.A. (gris claro) y topografía de la zona. 


Zona susceptible a inundarse 


Figura 12. Mapa Nacional Mexicano de Índice de Inundación. 
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posible identificar estos cuerpos de agua es 
relevante, pues la metodología no ha incorpo- 
rado dicha información. 

La figura 12 también muestra una amplia- 
ción parcial de la Región Hidrológica Grijalva- 
Usumacinta, donde se puede apreciar que 
el mapa nacional aparentemente indica que 
prácticamente toda la región está en riesgo de 
inundación debido al corto espacio para su 
impresión, aunque en realidad no ocurre así. 


Conclusiones 


Seha propuesto un índice deinundación basado 
en el Índice Topográfico, que adicionalmente 
considera aspectos climatológicos, edafológicos 
y de uso de suelo. El índice tiene como objetivo 
identificar zonas perennes y efímeramente 
saturadas o inundadas. El análisis realizado en 
la Región Hidrológica Lerma-Santiago muestra 
la efectividad del índice en la identificación de 
zonas perennemente inundadas, como lagos 
y ríos. Además, el índice guarda una relación 
razonable con el mapa de tasas de erosión y el 
modelo hidrológico-hidráulico de Agroasemex, 
S.A., en la planicie de inundación de Tabasco. 

El índice ha sido clasificado de manera 
binaria debido a su naturaleza bimodal. Se 
presenta un mapa nacional que indica zonas 
perennes o efímeramente inundadas y zonas 
libres de dicha amenaza. Los resultados son 
satisfactorios, aunque el tamaño de la impresión 
impide observar detalles. 

Cabe mencionar que el CEM de INEGI 
tiene severos problemas de continuidad. Los 
errores son prácticamente imperceptibles 
cuando se analizan a escala nacional, pero los 
análisis regionales vuelven evidentes dichas 


fallas. Los datos del SRTM están libres de las 
discontinuidades, y aunque su resolución es 
relativamente mayor, se considera apropiada 
para análisis de inundación de grandes 
extensiones, como el mapa nacional aquí 
presentado. 
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Abstract 


URIBE-ALCÁNTARA, E.M., MONTES-LEÓN, M.A.L. €: GARCÍA-CELIS, E. National 
Flood Index Map. Water Technology and Sciences, formerly Hydraulic engineering en 
Mexico (in Spanish). Vol. 1, no. 2, April-June, 2010, pp. 73-85. 


Floods are associated with factors such as: 1) river overflow; 2) flash floods, 3) high tides 
associated with hurricanes; and 4) breakage of control structures. Unfortunately, historical 
records of floods are practically inexistent. Therefore, risk assessments are generally based on 
indirect methods such as physical-mathematical modeling. The development of a National 
Flood Index Map for the identification of regions prone to floods is presented. The index 
proposed is based on the topographic index developed by Beven (1979) but additionally 
considers soil, hydrological and climatological factors. The process is detailed for the hydrologic 
region Lerma-Santiago. In addition, the results are compared with corresponding simulations 
of a routing model for the floodplain of Tabasco. The results indicate that the index is able 
to capture reasonably well the perennial and ephemeral flooded regions. The Mexican Flood 
Index Map based on the proposed methodology is included. 


Keywords: flood, index, risk, Mexico, Grijalva-Usumacinta, Lerma-Santiago. 
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